Ferromagnetic transition in a double-exchange system 
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We study ferromagnetic transition in three-dimensional double-exchange model. The influence of 
strong spin fluctuations on conduction electrons is described in coherent potential approximation. In 
the framework of thermodynamic approach we construct for the system "electrons (in a disordered 
spin configuration) + spins" the Landau functional, from the analysis of which critical temperature 
of ferromagnetic transition is calculated. 
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I. INTRODUCTION 

The recent rediscovery of colossal magnetoresistance 
(CMR) in doped Mn oxides with perovskite structure 
Ri_a;Da;Mn03 (R is a rare-earth metal and D is a diva- 
lent metal, typically Ba, Sr or Ca) substantially in- 
creased interest in the double-exchange (DE) model . 
Several approaches were used lately to study the thermo- 
dynamic properties of the DE model, including the Dy- 
namical Mean Field Approximation (DMFA) (see [Q and 
references therein), Schwinger bosons and Variational 
mean- field approach Calculated critical tempera- 

ture Tc of the ferromagnetic (FM) transition in three- 
dimensional DE model, surprisingly well agrees with the 
results of Monte Carlo method , and high tempera- 
ture expansion [|| for all values of the itinerant electrons 
concentration. Still we decided to present our calculation 
of the Tc- We believe that our derivation is interesting 
by itself; it also opens the opportunity for the analysis 
of the concurrent action of the impurity potential and 
magnetic disorder [p 10|. 



II. HAMILTONIAN AND CPA EQUATIONS 

We consider the DE model with the inclusion of the 
single-site impurity potential. We apply the quasiclassi- 
cal adiabatic approximation and consider each core spin 
as a static vector of fixed length 5 (Sj = S'n.i, where rij 
is a unit vector). The Hamiltonian of the model i?([ni]) 
in site representation is 



+ (e,: - Jrii • a) 5, 



(1) 



where ti-j is the electron hopping, is the on-site en- 
ergy, J is the effective exchange coupling between a core 
spin and a conduction electron and a is the vector of the 
Pauli matrices. The hat above the operator reminds that 
it is a 2 X 2 matrix in the spin space (we discard the hat 
when the operator is a scalar matrix in the spin space). 



The Hamiltonian (|^) is random due to randomness of a 
core spin configuration [rij]; we also take into account the 
randomness of the on-site energies e;. 

To treat Hamiltonian (h ) we use the coherent potential 



Approximation (CPA) i.e. we present it as 



H — I ti^j 



^5, 



(. 



Jn, • (7 - S 5, 



(2) 



(the site independent self-energy S(i?) is to be deter- 
mined later), and construct a perturbation theory with 
respect to random potential V . To do this let us intro- 
duce the T-matrix as the solution of the equation 



T = V + VGoT, 



(3) 



where Gq = {E — H^) ^ . For the exact Green function 
G we get 



G — Go + GqTGq 



(4) 



The coherent potential approximation (CPA) is ex- 
pressed by the equation 



G 



Go 



0, 



where Ti is the solution of the equation 
f, = V, + V,g{E-t)f,, 

and 

No{e) 



g{E) = (Gom^ 



E~e 



de, 



(5) 



(6) 



(7) 



where Nf){£) is the bare density of states. Equation for 
the T-matrix can be presented as an algebraic equation 
for the 2x2 matrix E: 



1 - V,g{E - S) 



0. 



(8) 



Eq. (0) coincides exactly with that of Ref. p 
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III. LANDAU FUNCTIONAL OF A 
DOUBLE-EXCHANGE SYSTEM 

In this part we start from the exact partition function 
of the system electrons + spins 



exp 



-/3i^e([n,]) + /3H^n, 



Dn 



3 ' 



(9) 



where N is the number of core spins, H is magnetic field, 
(3 = 1/T, and -Fe([nj]) is the free energy of the electron 
subsystem for a given core spin configuration [rij]. Our 
aim is to obtain Landau functional of the system. 

In the mean-field approximation all energy levels 
depend only upon macroscopic magnetization m = 
74 S j ■ (Mean-field means that we do not take into 
account large scale fluctuations of macroscopic magne- 
tization, but we certainly take into account microscopic 
fluctuations ofn^.) Hence 

Fe([n,]) « i^e(m), (10) 

and the partition function can be written as 

Z = y exp [-/3F(m, H)] dm, (11) 

where 

F{m.,U)=Fe{ra)-Nnm-TS{m)] (12) 
the quantity 



S{m) = In / 5 



TV ^ 

i 



(13) 



is the entropy of the core spin subsystem. We can identify 
the exponent with the Landau functional of the whole 
system Further on we'll be interested in the case 

H = 0; in this case F{m., H) = F(m). 

At high temperatures the minimum of the functional 
is at m = 0. At low temperatures the point in — Q cor- 
responds to the maximum of the functional. So if we 
expand F,, and S with respect to ni 



F^(m) ^ Fq- Fam^ 
S{m) = So — S2m^ 



(14) 



the critical temperature Tc is found from the equation 



F2 - T,S2. 



(15) 



In fact, we need only the quadratic terms of the Landau 
functional to find the Tc. 

Here it is appropriate to note that in the approach of 
Ref. ijl^ the thermodynamic part of the solution is based 
on the equation 



(n) = j dn^P{n)n, 



(16) 



where P{n) is the Boltzmann factor for a spin, deter- 
mined from the DMFA ansatz. 

To calculate the entropy we expand the (5-function in 
Eq. (|l^) in Fourier integral 



(17) 



which gives an opportunity to calculate integrals with re- 
spect to Duj . After simple algebra we obtain (ignoring 
irrelevant multipliers in the argument of the logarithm) 

N ^ 



S{m) = In <^ (47r) 



N 



sm{Nrm) 



rdr 



(18) 



Calculated integral by the steepest descent method, up 
to the terms of order of we obtain 



S = N 



ln(47r) 



3m 



(19) 



The grand canonical potential for the electron subsys- 
tem is 

f + 00 



N 

— I3n{ii, m) = — 

TT 



In 



1 + e"'^(^"'') 



liJi gc{E)dE, 
(20) 



where gc is the charge locator gc — Tt g. The free energy 



Fe{N, m) = m) + fiN; 



(21) 



the connection between the number of holes per cite x 
and chemical potential /i is given by the equation 



l-. = i 

TT 



fiE)lm gc{E)dE, 



(22) 



where f{E) is the Fermi distribution fmiction. 

The calculation is especially simple, when we can con- 
sider the electron gas as degenerate and ignore the terms 
of order of T/zi (that is we can consider the temperature 
of electron gas as being effectively equal to zero). In this 
case 



N f 

Fe = — / E Imgc{E)dE. 



Expanding gc with respect to m 



3c .90 + 92m + 



we obtain 



F2 = - (po- E)lmg2iE)dE, 

J-co 

where /iq is found from the equation 
1 



1 - a; = - 

TT 



Im go{E)dE. 



(23) 
(24) 
(25) 

(26) 
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IV. Tc FOR SEMI-CIRCULAR DENSITY OF 
STATES 

We consider semi-circular (SC) bare (i.e. for e = and 
J = 0) DOS 



at |e| < W and A^oC^) = otherwise, for which 
9{E) 



2 

W 



E 
W 



(27) 



(28) 



Hence we obtain 

t^E-2wg-g-\ (29) 
where w = Thus, Eq. (^ can be presented as 

' (E - e + Jna ~ 2wg)'^) . (30) 



(31) 



It is convenient to write the locator g in the form 
9 = ^(ffcZ + gsO-), 



where / is a unity matrix. For the charge locator gc and 
spin locator we obtain the system of equations 



{E - e- wgc) - ( Jn + wg^ 



Jn + w^s 



[E - t - wgc) - ( Jn + wgs 



(32) 



In the strong Hund coupling limit ( J — s- oo) we obtain 
from Eqs. (32) two decoupled spin sub-bands. For each 
sub-band, after shifting the energy by ± J we obtain 



9c 



9s 



E — e — wgc ~ wn^g_ 
E — e — wgc — wn^gs 



(33) 



where gs = (0,0, ^s) (axis OZ is directed along m). 

For a saturated FM phase {m — 1), we obtain = gc, 
and closed equation for the charge locator is 



9c 



1 



E -e- 2wg 



(34) 



For a paramagnetic (PM) phase (m = 0), we obtain 
gi = 0, and closed equation for the charge locator 
gdm = 0) = go is 



90 = 



1 



E — e — wgo 



(35) 



In both cases only the averaging with respect to random 
on-sight energies is left. So, the DOS in PM phase is 
equal to the DOS in FM phase of another model, with 
decreased by a factor of bare bandwidth. 

From here on we ignore the on-site energies disorder 
and put e = 0. In this case the solution of Eq. ( ^5| ) for 
go coincides with Eg . d2q ), with W/V2 substituting for 
W. Expanding Eq. (|33[) with respect to m we obtain the 
system 



,w- 



92m = goiwg2 + mwgs)+go^9s^ 

2 ^ 

gs = mgo + go-^9s- 



(36) 



Substituting the value of 52 into Eq. (g^) , we obtain the 
critical temperature: 



Tc 



2w 
3^ 



(^0 - -E)Im 



9l 



a-Wo){i-f9iy 



dE. 



(37) 



Calculated integral analytically, we get for the Tc'. 



Tc 



V2 
An 



W (38) 



where the reduced chemical potential y 
found from the equation 



1 1 

2^ n 



sin ^ y + y\Jl-y'^ 



^/2^o/W^ is 



(39) 




0.25 0.5 0.75 

X (holes per site) 



Fig. 1. The calculated Curie temperature for the semi-circular 
DOS at J = 00 as function of hole concentration x. 

One notices the particle-hole symmetry, namely 
Tc{x) = Tc{l -x). At X ^ and x ^ 1, 0. The 

optimal doping turns out to be 0.5, with Tc = 0.045M^. 
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V. DISCUSSION 

Any calculation of the Tc in the DE model includes, on 
one hand, calculation of the energy of electron interact- 
ing with randomly directed core spins, and, on the other 
hand, the analysis of phase transition in the system of 
spins with the interaction inferred from the solution of 
the electron part of the problem. The difference between 
our approach and that of Rcf. , is similar to the differ- 
ence between the two mean-field approximations for the 
Heisenberg model: one based on the analysis of a local 
spin in a coherent field, which by itself depends upon the 
averaged magnetization, the other based on the expan- 
sion of the partition function . For the Heisenberg 
model the two approaches give exactly the same result. 
This turned out to be the case also the case for the DE 
model. Our result for the Tc (Eq. (|37|)), using integra- 
tion by parts, can be exactly reduced to Eq. (49) of Ref. 
||14|, which in our notation is 



2w 


/ Im 


[ d 1 




' — OO 





dE. 



(40) 



The result of Refs. |^ for the Tc in our notation is 

Tc=^{l-y'f'w. (41) 



It is close to our result numerically, but the functional 
dependence of the Tc upon hole concentration is essen- 
tially different. The difference comes from our use of the 
CPA for the description of electrons, in distinction from 
the mean-field decoupling of charge and core spin fluctu- 
ations in de Gennes' spirit, used in Refs. ||]. (In the part 
of Refs. [^1 where the charge-spin coupling is taken into 
account, the explicit result for the Tc is not presented.) 

The DE scenario of the FM transition in manganites 
attracted especial attention after its validity was ques- 
tioned in Ref. [|l6|. In the following discussion, the posi- 
tion of one side can be expressed by reformulation of the 
rhetorical question by Hubbard : " How, in the itiner- 
ant model, explain a Curie temperature ^ 1000° for iron, 
when calculations always give an exchange field ~ 1 — 2 
eV?" in the form: How, in the DE model, explain a Curie 
temperature 300° for a manganite (at optimal doping), 
when calculations always give a bandwidth ~ 1 — 2 eV? 
The fallacy is based on the implicit assumption that the 
Tc should be of the order of the energy difference be- 
tween the fully ordered FM state and the paramagnetic 
state (which is the smallest of the band width and the 
exchange integral). It does not take into account that the 



fluctuations which restore magnetic symmetry in itiner- 
ant models are local fluctuations in the orientation of 
magnetic moments, with much lower energy 

[0- That 

is why the numerical coefficient in Eq. ( |3^ ) turns out to 
be much less than one. 
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